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Abstract. In this article, we review the prospects for the Fermi satellite (formerly 
known as GLAST) to detect gamma rays from dark matter annihilations in the Cen- 
tral Region of the Milky Way, in particular on the light of the recent astrophysical 
observations and discoveries of Imaging Atmospheric Chcrenkov Telescopes. While 
the existence of significant backgrounds in this part of the sky limits Fermi's discovery 
potential to some degree, this can be mitigated by exploiting the peculiar energy spec- 
trum and angular distribution of the dark matter annihilation signal relative to those 
of astrophysical backgrounds. 
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1. Introduction 

Despite the numerous cosmological and astrophysical indications of the presence of non- 
baryonic dark matter (DM), the particle nature of this substance remains unknown. 
If the DM consists of weakly interacting massive particles (WIMPs), an important 
tool for inferring their properties could be "indirect detection", i.e. astrophysical 
observations of the annihilation (or decay) products of DM in our Galaxy or beyond. 
For WIMPs with masses at or around the electrowcak scale, mx ~ 0(0.1-1) TeV, the 
annihilation products are typically found at GeV-TeV energies, the domain of high- 
energy astrophysics. Of the different annihilation products, gamma rays and neutrinos 
have the important advantage of retaining directional information while not suffering 
energy losses. The very small cross sections of neutrinos, however, make their flux 
from the region of the Galactic Center very difficult to detect. On the contrary, the 
gamma-ray spectrum from DM annihilation or decay may be detectable with sufficient 
statistics, energy resolution, and over an extended angular distribution, to provide a 
very distinctive set of information related to both the particle identity of the WIMP 
and its astrophysical distribution. A major challenge is in separating the DM signal 
from any astrophysical backgrounds, whose energy spectrum and angular distribution 
are not well known. 

A major change in the prospects for DM detection has occurred in the last few 
years, following the discovery of a bright astrophysical source of TeV gamma rays from 
the Galactic Center. As a result of this discovery, we now know that DM emission from 
the Galactic Center will not be detectable in a (quasi) background-free regime, and — 
unless one turns the attention to other targets — the peculiar spectral shape and angular 
distribution of the signal must be used to extract it from this and other backgrounds. 
This topic is the main subject of this review. The remainder of this paper is structured 
as follows: In Sec. 2 we briefly describe the capabilities of current high energy gamma- 
ray telescopes. In Sec. 3 we review the features of the DM annihilation signal, both in its 
energy (Sec. 3.1) and angular dependence (Sec. 3.2). In Sec. 4 we review the prospects 
for detecting or constraining DM annihilations in the Galactic Center region, while in 
Sec. 5 we briefly discuss some interesting aspects of the diffuse signal from the inner 
halo. In Sec. 6 we report our conclusions. Appendix A briefly illustrates the qualitative 
changes to the considerations we develop in the main text when decaying (as opposed to 
annihilating) DM is considered. Appendix B describes a subtlety regarding the angular 
shape of the signal arising for a velocity-dependent annihilation cross section. 

2. The Instruments 

The ability of gamma-ray experiments to identify DM annihilation radiation from the 
Galactic Center region relies on the effective area, angular and energy resolution of 
the existing telescopes, as well as the rejection of other (mostly hadronic) cosmic ray 
events contaminating the gamma-ray sample. Since the mean free path of gamma rays 
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is much shorter than the atmospheric slant depth, direct observations in the GeV region 
and above can only be done from space — which is the strategy pursued by the LAT 
detector on the Fermi gamma-ray space telescope (formerly GLAST) [1] — or indirectly 
by ground-based Imaging Atmospheric Cercnkov Telescopes (lACTs) such as HESS [2], 
MAGIC [3], VERITAS [4] and CANGAROO-III [5]. In the latter category, the direction 
and energy of the primary particle hitting the atmosphere is reconstructed from the 
Cherenkov emission of the secondary charged particles generated in the atmospheric 
shower. 

These differences lead the two classes of experiments to adopt different strategies 
in the search for DM. Fermi-LAT is very effective in rejecting hadronic events, and 
continuously monitors a large fraction of the sky, but has an effective area of only 
~ Im^, far smaller than that of ground-based telescopes, ~ 10^ m^. On the other 
hand, lACTs study small angular fields and have a lower rejection capability, but much 
greater overall exposure. As a consequence, diffuse gamma-ray signals are better probed 
by Fermi-LAT. Any unidentified sources detected by Fermi-LAT which lack a low-energy 
counterpart could be potentially attributed to DM substructure. lACTs would be very 
effective in providing detailed follow-up observations of such sources. 

In addition, the accessible energy range is very different between these two classes 
of experiments: ~100 MeV to 300 GeV for Fermi-LAT, and above ~100 GeV for ACTs. 
This difference makes Fermi-LAT most sensitive to DM particles lighter than a few 
hundreds GeV, while lACTs are better suited for TeV-scale or heavier WIMPs. On 
the other hand, the Fermi-LAT has poorer angular resolution than lACTs, so it is 
less accurate in the localization of point-like sources. For both instrument classes, the 
search for indirect DM signatures is among the top physics priorities. The reach of 
Fermi-LAT and current and future lACTs has been recently assessed in Ref. [6] and 
Ref. [7], respectively. 

3. The Dark Matter Signal 

The differential flux of gamma rays (photons per unit area, time, energy and steradian) 
produced in DM annihilations! is described by 



where {av) is the WIMP annihilation cross section multiplied by the relative velocity of 
the two WIMPs (averaged over the WIMP velocity distribution), mx is the mass of the 
WIMP, p is the position-dependent DM density, and the integral is performed over the 
line-of-sight (los) in the direction of the sky, ft. The gamma-ray spectrum generated per 
WIMP annihilation is dN^/dEj, it has units of Energy"^ and its integral over energy 
is equal to 1. If the DM is not its own antiparticle as assumed here, Eq. (1) should 
be multiplied further by a factor 1/2 (if X and X are equally abundant). The factor 

I See Appendix A for the case of DM decay. 




(1) 
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in square brackets in Eq. (1) depends only on particle physics: in particular, cross 
section, mass, and the spectrum of gamma rays produced through DM annihilations 
depends on the nature of the WIMP. The integral over the line-of-sight determines 
instead the angular dependence of the signal and is controlled by the astrophysical 
distribution of DM. We shall discuss each of these two terms in the following, assuming 
that this factorization holds (see Appendix B for a discussion of possible violations of 
this hypothesis). 

In convenient units, Eq. (1) can be recast as: 

cm-2 s-i sr-i ^ ^ dEj^ -ph \ mx J 

The dimensionless function J{Q) depends only on the DM distribution in the halo and 
is defined by convention as [8] 

For typical halo models (see Sec. 3.2) this is a function strongly peaked towards the 
Galactic Center (for an illustration, see for example Fig. Al). 

The benchmark value for the cross-section, (av) ~1 picobarn (or in cgs units 
~ 3 X 10^^*^ cm'^/s), is motivated by the fact that a WIMP annihilating with such a 
cross section during the freeze-out epoch will be generated as a thermal relic with a 
density similar to the measured DM abundance (for a review, see Ref. [9]). WIMPs 
constituting the cold DM annihilate in the non-relativistic limit. If annihilations take 
place largely through S-wme processes, then the annihilation cross section of WIMPs 
in the Galactic halo (i.e. in the low velocity limit) will also be approximately equal to 
this value, which justifies the benchmark value used in Eq. (2). Yet, it is important 
to stress that much lower signals are possible (e.g. if P— wave annihilation dominates 
the freeze-out process), as well as significantly enhanced ones in models where the DM 
is non-thermally produced in the early universe (see e.g. Ref. [10]). Some further 
considerations and refs. can be found in Appendix B. 

It is also worth noting that for S-wave annihilating thermal relics the indirect 
detection signal has two advantages compared to direct detection via nuclear recoils 
in underground detectors: i) it is proportional to a relatively large annihilation cross 
section; ii) it is less dependent from the particle physics details. For example, even 
under the (possibly unrealistic) assumption that DM annihilates mostly into quarks of 
the first generation, the natural expectation value for DM-nucleon elastic cross section 
is at the level of {mN/mxY x Ipb^ 10"''° (mx/0.1 TeV)^ cm^. For DM annihilating 
predominantly into heavier particles, a further suppression is expected. 

3.1. Energy Spectra 

The gamma-ray spectrum from DM annihilation originates from several different 
contributions. Typically, the most abundant source of photons is the hadronization 
and/or decays of unstable particles. For example, neutralinos in the Minimal 
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Figure 1. The gamma-ray spectrum per annihilation for a 100 GeV (left) and 500 
GeV (right) WIMP. Each curve denotes the result for a different dominant annihilation 
mode. From Ref. [12]. 



Supersymmetric Standard Model (MSSM), dominantly annihilate to final states 
consisting of heavy fermions bb, tt, t^t~ (i.e. bottom quarks, top quarks and tau 
leptons, respectively) or bosons ZZ, W~^W~, HA, HA, ZH, Zh, ZA, W^H"^, where 
W^, Z are the gauge bosons mediating the weak interactions and H, h, A and are 
the Higgs bosons of the MSSM [11]. With the exception of the r+r^ channel, each of 
these annihilation modes result in a very similar spectrum of gamma rays, dominated 
eventually by the decay of mesons, especially 7r°, generated in the cascade. In Fig. 1 we 
show the predicted gamma-ray spectrum per annihilation, for several possible WIMP 
annihilation modes. The harder spectrum here shown from annihilation into the r+r~ 
channel is not typical of SUSY models, although it might be a distinctive feature of DM 
annihilating dominantly to charged leptons pairs; this arises, for example, in Kaluza- 
Klein DM models (for a review of particle DM models, see e.g. [13]). 

Although these secondary photons provide the dominant emission, other important 
channels can exist at the one-loop level. A particularly striking signature would be the 
mono-energetic photons resulting from final states such as 77, 7Z or 7/1 (see Ref. [14] 
for a discussion of these processes within the context of supersymmetry). Unfortunately, 
such processes are expected to produce far fewer events than continuum emission and 
in typical models can not easily be detected (for a counterexample, see Ref. [15]). 
Bremsstrahlung (with an additional photon appearing in the final state) is automatically 
present when annihilations produce charged final states, and can dominate the high 
energy region of the spectrum when those charged annihilation products are much lighter 
than the WIMPs. This is particularly important when the tree-level processes to a pair 
of light fermions are disfavored by "selections rules" , but no suppression is present for 
three body final states [16, 17]. For more details, see the review by L. Bergstrom in this 
issue [13]. Note that gamma-ray emission from DM annihilating into standard model 
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particles is unavoidable, even in the most extreme case when only neutrino final states 
are allowed at tree level, due to W and Z-strahlung [18, 19]. 

3.2. Angular Shape Of The Signal 

The gamma-ray spectrum and angular distribution predicted by Eq. (1) is rather 
general (sec, however. Appendix B), and could be applied to the case of a smoothly 
distributed Galactic halo, or alternative targets such as dwarf spheriodal galaxies, 
microhalos/clumps, density spikes around intermediate mass black holes, or the 
integrated extragalactic diffuse background (described in other contributions to this 
special issue). While some of these targets have interesting observational prospects, 
the intensity of these signals is strongly dependent on often unknown cosmological and 
astrophysical properties, such as the quantity of small-scale structures in DM halos or 
the population of intermediate mass black holes in the Milky Way. In this article, we 
will focus on the smooth Galactic halo. 

Naively, by inputing the DM density profile inferred from kinematical observations 
of the Milky Way into Eq. (1), one could obtain an approximate lower bound on the 
gamma-ray flux from DM annihilation. In principle, this would enable us to translate the 
observations from ground or space-based gamma-ray telescopes into constraints on the 
particle physics properties of the WIMP (mass, annihilation cross section and dominant 
modes). Unfortunately, even the average, smooth distribution of the DM particles in 
our Galaxy is not well known, especially in the volume within the solar circle. This is 
due to the fact that the baryonic material dominates the gravitational potential in the 
inner Galaxy, and lacking a detailed knowledge of its distribution, a reconstruction of 
the DM distribution "by subtraction" is unfeasible. It is not surprising, then, that very 
different profiles have been claimed to fit the observations (for example, compare the 
results of Rcf. [20] and Ref. [21]). 

On general grounds, a class of spherically symmetric, smooth halo distributions 
can be used to approximately fit both the observed rotation curves of galaxies and the 
results of numerical simulations of DM halos. A non-trivial angular dependence of the 
gamma-ray signal results from the off-center position of the Sun within the halo (see e.g. 
Ref. [22] and references therein for a discussion of sub-leading effects determining the 
angular distribution of the signal). The function J introduced in Eq. (2) and Eq. (3) then 
depends only on the angle 6 between the observed direction of the sky and the Galactic 
Center or, in terms of galactic latitude b and longitude /, only on cos^^ = cos 6 cos/. The 
radial variable, r, can be expressed in terms of the relevant quantities {i, 6} as 



where r© fa 8.33 ± 0.35 kpc [23] is the distance of the Solar System from the Galactic 
Center. Typically, one considers a distribution of the form 




(4) 




(5) 
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where po is a normalization constant and is a characteristic radius below which the 
profile scales as r"'''. A very well known, universal profile of this class fit to DM-only 
(i.e. neglecting baryons) N-body simulations has been proposed by Navarro, Frenk and 
White (NFW) [24], corresponding to the choice {«,/?, 7} = {1.0,3,1.0}. Steeper or 
softer profiles have also been extensively discussed, such as that proposed by Moore 
et al. [25] and Kravtsov at al. [26], respectively. While simulations (and data as well) 
typically agree on the shape of profile in the outskirts of the halos, a disagreement clearly 
exists concerning the inner slope, 7. More recent simulations [27, 28, 29, 30, 31] suggest 
that halo density profiles are better represented by a function with a continuously- 
varying slope, as the one proposed by Einasto [32] 

p(r) = p_2e "LV-^; J, (6) 

with clear hints for non-universality in the form of halo-to-halo variations in the quantity 
a. The quantity r_2 — 25 kpc denotes the radius at which the logarithmic slope of the 
profile, dlogp/dlogr, assumes the value —2; the other free parameters are a and the 
overall normalization, here chosen as the density p_2 = p(r_2). At r ^ 1 kpc, these 
newly proposed fitting formulae provide only marginal improvement with respect to the 
more traditional ones. At smaller radii, however, these recent results lead us to expect 
the inner slope of DM halos to be shallower that that predicted by Moore et al., and 
probably shallower than NFW, as well. 



Model 


a 




7 


Pe 

[GeVcm-3] 


[kpc] 


Moore et al. 


1.5 


3 


1.5 


0.27 


28 


NFW 


1.0 


3 


1.0 


0.30 


20 


Kravtsov 


2.0 


3 


0.4 


0.37 


10 



Table 1. Parameters describing some common halo profiles of the form described by 
Eq. (5), where is the DM density at the Solar distance from the GC. See text for 
more details. 



There are a number of astrophysical processes that may potentially modify the DM 
distribution, none of which are taken into account in the above-mentioned profiles. It is 
very difficult to reliably account for these effects in simulations, and only a few results 
are available (see e.g. Refs. [33, 34, 35]). Qualitatively, since baryons can cool and 
contract, one expects them to steepen the gravitational potential in the central regions 
galaxies and, as a result, enhance the DM density [36]. On the other hand, other 
feedback or frictional effects have been proposed that could reduce the DM density in 
the inner halos and bring the prediction closer to observations (for a recent review, see 
Ref. [37]). For some galaxies, it has actually been argued that flat-cored profiles fit the 
observations better than cusped profiles. 

Even greater uncertainties exist concerning the DM distribution at the very center 
of the Milky Way, in the region immediately surrounding the central supermassive black 
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hole. Adiabatic accretion may lead to the formation of a spike in the DM distribution, 
resulting in a very high DM annihilation rate in the innermost parsecs of the galaxy [38] . 
Mergers as well as scattering on the dense stellar cusp around the central black hole may 
potentially destroy density enhancements, however. In general, these effects only affect 
regions too close to the Galactic Center to be resolved angularly by present detectors, 
leaving only the energy spectrum to be used for separating the DM signal from the 
background. Yet, such a spike might lead to a measurable gamma-ray flux from the 
innermost angular bin, even in presence of relatively large astrophysical backgrounds. 
For more details, we direct the reader to the references cited in Ref . [9] or Ref . [39] . 

Given these considerably uncertainties, wc have chosen to use three illustrative 
choices for the DM halo profiles: the Moore et ai, NFW, and Kravtsov profiles (with 
the parameters given in Table 1). In each case, as in Ref. [40], the normalization has 
been chosen so that the mass contained within the solar circle provides the appropriate 
DM contribution to the local rotational curves. In all of the cases we have discussed, 
even for the most conservative cored profile, the DM signal peaks at the Galactic Center. 
In presence of an isotropic (or even vanishing) astrophysical background, the Galactic 
Center region thus becomes the natural location to look for a DM signal. Unfortunately, 
there is also a higher background density as we look toward the inner region of our 
Galaxy. Still, as a first step, one can consider the detection prospects for this region, 
a task wc address in Sec. 4, before turning to more general arguments regarding the 
optimal regions of the sky for DM searches in Sec. 5. 

4. The Galactic Center 

The Galactic Center has long been considered to be among the most promising targets 
for the detection of DM annihilation, particularly if the halo profile of the Milky 
Way is cuspcd in its inner volume [8, 41]. This has been complicated, however, by 
the recent discovery of astrophysical gamma-ray sources from the Galactic Center. 
Following an earlier claim by the WHIPPLE lACT [42], very high-energy gamma 
rays from the Galactic Center have been detected by HESS [43], MAGIC [44], and 
CANGAROO-II [45]. This source is consistent with point-like emission and is located 
at I = 359°56'41.1"±6.4" (stat), b = -0°2'39.2"±5.9" (stat) with a systematic pointing 
error of 28" [46], coincident with the position of Sgr A*, the black hole constituting the 
dynamical center of the Milky Way. The spectrum of this source is well described by 
a power-law with a spectral index of a = 2.25 ± 0.04(stat) ± O.lO(syst) over the range 
of approximately 160 GeV to 20 TeV. Although speculations were initially made that 
this source could be the product of annihilations of very heavy (~10 to 50 TeV) DM 
particles [47], this interpretation is disfavored by the power-law form of the observed 
spectrum and the wide energy range over which it extends (see Fig. 2). The source 
of these gamma rays is, instead, likely an astrophysical accelerator associated with our 
Galaxy's central supermassive black hole [49]. In recent analyses, this source has been 
treated as a background for DM searches [50, 51] . Given the presence of this background. 
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Figure 2. Spectral energy density of ttie Galactic Center source as measured by HESS 
in 2004 (full points) and 2003 [43] (open points). Upper limits are 95% CL. The shaded 
area shows the power-law fit dN/dE - E-^ , with T = 2.25 ± 0.04 (stat.)±0.10 (syst.). 
The dashed line illustrates typical spectra of phenomenological MSSM DM annihilation 
for best fit neutralino masses of 14 TcV. The dotted line shows the distribution 
predicted for Kaluza-Klein DM with a mass of 5 TeV. The solid line gives the spectrum 
of a 10 TcV DM particle annihilating into r+r" (30%) and bb (70%). From Ref. [48]. 



the prospects for detecting DM annihilation products from the Galactic Center appear 
considerably less promising than they had a few years ago. 

The main challenge involved in DM searches with Fermi will be to distinguish 
the signal from this and other backgrounds. This is made particularly difficult by 
our ignorance regarding the nature of these backgrounds. The Galactic Center is 
indeed a complex region of the sky at all wavelengths, the gamma-ray window being no 
exception [52]. 

An attempt to use angular and spectral information to separate DM annihilation 
products from these backgrounds was performed in Ref. [12], studying a 2° x 2° 
region arond the Galactic Center (see also the analysis of the Fermi DM team in 
Ref. [6]). This study considered two known point-source backgrounds: a yet unidentified 
source detected by EGRET approximately 0.2° away from the dynamical center of our 
galaxy [53, 54], and the lACT source discussed above. In addtion, a diffuse spectrum 
with a free power-law index was also included. More in detail, we describe the spectrum 
of the source revealed by lACTs at the Galactic Center as a power-law given by: 

. P s -2.25 

= 1.0 X 10-^ (■^] GeV-i cm-2 s'^. (7) 
The flux from the EGRET source slightly off-set is instead modeled as: 

<l>^^ = 2.2 X 10-' ( GeV-^ cm'^ s"^ , (8) 
VGeV/ 

where the exact value of the cutoff energy (here 30 GeV) is somewhat arbitrary, but 
reflects the fact that this source has not been observed yet by lACTs a.t E > 100 GeV. 
Finally, we allow for a diffuse/unresolved flux with spectrum 

$'^'«^(Aa) = ^(^) °GeV-^cm-2s-^r-\ (9) 
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Figure 3. The projected exclusion limits at 95% confidence level from Fermi-LAT 
(after ten years) on the WIMP annihilation cross section, as a function of the WIMP 
mass. The region above the dotted line is already excluded by EGRET [53]. The solid 
and dashed lines show the projections for Fermi-LAT for an assumed isotropic diffuse 
background and the limiting case where the astrophysical background has exactly the 
same angular distribution of the DM signal, respectively. In the left and right frames, 
the NEW and Moore et al. halo profiles have been adopted, respectively. Also shown 
are points representing a random scan of supersymmetric models. From Ref. [12]. 



where a is allowed to vary between 1.5 and 3.0. We adopt an overall normalization, A, 
such that the integrated flux of the diffuse background between 1 GeV and 300 GeV in 
a 2° X 2° field of view around the Galactic Center is equal to 10~^ cm~^ s~^ sr~^. We 
do not, however, assume that this normalization is known in our analysis, leaving open 
the possibility that some of the diffuse gamma rays observed are the product of dark 
matter annihilations. 

A multi-parameter analysis of the simulated sky against models including a 
contribution from DM annihilation radiation yields the projected exclusion limits at the 
95% confidence level shown in Fig. 3 (for ten years of collection time by Fermi-LAT). 
The simulated sky used to produce this figure contains the two resolved sources and the 
diffuse background but no dark matter, and to derive the exclusion limits we compare 
the signal thus obtained to a model which includes both the backgrounds and a signal 
from dark matter. The left and right panels refer to the cases of an NFW profile and 
Moore et al. profile (replaced by a fiat core within 10~^ pc of the Galactic Center to avoid 
a divergency), respectively. The WIMP is always assumed to annihilate dominantly to 
bb. This assumption is not particularly restrictive, since for many annihilation modes 
the spectrum would look very similar (see Fig. 1). Apart for some cases, as possibly 
Kaluza Klein DM, the gamma-ray spectrum alone does not help much in discriminating 
among several DM candidates, see e.g. [55] for details. 

To give a feeling for the numbers involved, the integral of the function J over the 
considered region leads to a factor ~ 1.3 for the NFW model; in this model, for the 
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Figure 4. The ability of Fermi-LAT to measure the annihilation cross section and mass 
of dark matter after ten years of observation. A benchmark scenario with mx = 100 
GeV, {(jv) = 3 X 10^^^ em'^/s and an NFW halo profile has been used. The inner 
and outer contours in each frame represent the 2 and 3 a regions, respectively. In the 
left frame, the halo profile shape was treated as if it is known in advance. In the 
right frame, a marginalization over the inner slope of the profile was performed. From 
Ref. [121. 



benchmark 100 GeV WIMP with (erf) ~ 3 x 10^^^ cm^/s one would expect less than 
3000 events from DM above 1 GeV, to be compared with more than 2.3 x 10^ background 
photons. The background is almost evenly split between the EGRET point source and 
the diffuse flux; the lACT source contributes less than 7000 events, but it is located 
at the GC and, differently from the EGRET source, is not cut-off at high energy, so 
it is important to include it especially for high WIMP masses. The solid line in each 
frame represents the limit found if the diffuse background is assumed to be distributed 
isotropically, while the dashed line represents the conservative limit obtained if the 
diffuse background has the same angular distribution as the DM signal {i.e. the case 
in which angular information is not useful in disentangling the signal from the diffuse 
background). For values of (av) below the corresponding lines, a pure background model 
is expected to be consistent with the data. The fact that the limits are significantly 
stronger in the uniform background case is the manifestation of the improved sensitivity 
which can be achieved by an analysis including both energy and angular information. 
For comparison, in Fig. 3 it is also shown the region already excluded by EGRET [53] 
(above the dotted line) and the mass and cross section of neutralino models found in 
a random scan over supersymmetric parameters, as calculated using DarkSUSY [56]. 
As expected, many of the models cluster around (av) ~ 3 x 10~^^ cm^/s, the value 
required of a thermal relic annihilating via an S-wave amplitude. Each point shown 
represents a model which respects all direct collider constraint and generates a thermal 
DM abundance consistent with the observed DM density. In the scan the SUSY 
parameters varied were M2, |/i| and up to 2 TeV, and up to 1 TeV and 
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Figure 5. The ability of Fermi-LAT to measure the inner slope of the halo profile and 
the mass of dark matter particle (marginalizing over the annihilation cross section) 
after ten years of observation. Here, a benchmark scenario with mx — 100 GeV, 
(av) = 3 X 10~^^ cm'^/s and an NFW halo profile has been used. The inner and outer 
contours represent the 2 and 3(7 regions, respectively. From Ref. [12]. 



tan/3 up to 60. Also, the gaugino masses were assumed to evolve to a single unified 
scale, such that Mi ^ O.SMa, Mg ^ 2.7M2. 

Should gamma rays be identified as having been produced in DM annihilations, 
such observations could then be used to measure the characteristics of the DM 
particle, including its mass, annihilation cross section and spatial distribution. Such 
determinations are an important step toward identifying the particle nature of DM. A 
calculation similar to the one leading to the results of Fig. 3 can be performed, this 
time including in the template a contribution from a fiducial model of DM annihilation, 
and asking the accuracy by which a reconstruction of its input properties is possible. In 
Fig. 4, the ability of Fermi-LAT to determine the WIMP mass and annihilation cross 
section for WlMPs distributed with an NFW halo profile is illustrated, for the case of 
an annihilation cross section of 3 x 10~^^cm^/s and a mass of 100 GeV. In each frame 
the projected 2 and 3 a constraints on the input parameters are reported, assuming 
an isotropic diffuse background (in addition to background point sources). In the left 
frame, we treat the shape of the halo profile (NFW) as if it is known in advance. Of 
course, this is not a realistic assumption, and a less accurate determination of the WIMP 
mass must be expected in a more realistic treatment. In the right panel of Fig. 4 we 
report the results obtained by marginalizing over the inner slope of the halo profile, 7. 
In the absence of an assumption on the inner halo slope, the constraint on the DM mass 
worsens by a factor ~2. 

If the spectrum and angular distribution of gamma rays from DM annihilations 
in the Galactic Center region are sufficiently well measured, it will also be possible to 
measure the underlying DM distribution. Fig. 5 shows the results from the right frame 
of Fig. 4 in the {mx, 7} plane, marginalized over the annihilation cross section. In 
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the above benchmark model, the inner slope of the halo profile can be determined at 
approximately the ~ 10% level. 

A few caveats are in order: We are possibly oversimplifying the spectral shape 
of the background, since we are extrapolating the known point-like source properties 
from lower and higher energies. Also, the bounds are slightly optimistic, in the sense 
that we are considering the longest plausible lifetime of Fermi (ten years) and that 
further effects degrading the angular resolution or accounting for dead-time may loosen 
the constraints by up to a factor~ 0{2). Finally, we have included only the HESS 
and EGRET sources (in addition to the diffuse background) in our analysis, but it is 
likely that other astrophysical point sources with different spectra will be discovered by 
Fermi-LAT. The above estimates should be quite realistic as long as a limited number of 
discrete sources will "contaminate" the inner Galactic Center angular bins, so that they 
could be removed effectively by the Fermi team. On the other hand, it is worth taking a 
lesson from the "sudden" discovery of an astrophysical source just at the Galactic Center: 
In the worst case, the Galactic center region might so crowded with (yet unknown) 
sources that the separation of background and signal might be degraded with respect to 
expectations based on present knowledge. Still, it is worth mentioning that the bulk of 
the statistical significance of the dark matter annihilation signal typically does not come 
from the inner 0.1° around the Galactic Center (where the lACT source dominates), 
rather from the surrounding angular region. It is then interesting to see what are the 
perspectives to detect a DM signal from a more extended region, which we address in 
the following. 

5. The Inner Halo 

The emission of radiation per unit solid angle from DM annihilation is expected to be 
maximized at the Galactic Center. Yet, geometric factors and the presence of point- 
like and diffuse backgrounds make the choice of the optimal window size a non-trivial 
problem [6, 7, 8, 57, 58, 59]. Unless the DM halo is very cuspy towards the Galactic 
Center (say, cuspicr than NFW), the optimal strategy for DM searches is never to 
focus on the inner sub-degree around the Galactic Center. Rather, a window size up 
to ~ 50° or more is preferred. The optimal shape of the window depends on the 
angular distribution of the signal and backgrounds, but also on the details of the analysis 
(like energy cuts, astrophysical foreground removal, etc.). It could be optimized from 
morphological studies of the low-energy emission measured by Fermi, but a circular 
annulus around the Galactic Center or a "rectangular" window in Galactic coordinates — 
with an inner rectangular mask — generally work fairly well. Here we illustrate this 
issue within a simplified model: Based on EGRET data, in the diffuse background one 
can identify an isotropic, extragalactic component oc E~'^'^ and a Galactic component 
scaling as ~ which is however dominant at GeV energies and peaks towards 

the Galactic Plane (see [8, 59] for details). The energy spectrum of the DM signal 
thus mostly enters the game in determining which one of the two spectra (of different 
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angular shape) dominates the background: the harder the DM spectrum, the closer the 
background is to an isotropic emission, the smaller the optimal angular window. Still, 
even in the quite extreme§ case of 70 GeV monochromatic lines, the optimal angular 
window is very extended. This is illustrated in Fig. 6 , where we plot the relative signal- 
to-noise [S/N) as a function of the maximum galactic latitude, 6max, for a Galactic 
region 0.4° < |6| < fomax, 0° < |/| < /max = ^max- The solid lines refer to an hypothetical 
gamma-line emission for the benchmark "model I" in the inert Higgs doublet scenario 
(see Ref. [15] for details) for three different halo profiles. The dot-dashed line refers to 
the continuum spectrum of a 100 GeV neutralino annihilating into W^W~ with an NFW 
distribution, where indeed we see that the optimal window is even more extended than 
for the line signal: this is due to the fact that the background at the typical energies 
of the continuum DM photons, E < mx/10, is mostly dominated by the Galactic 
unresolved background, which peaks towards the Galactic Plane/Center. To give a 




b„ax (deg) 

Figure 6. The relative signal-to-noise as a function of the window size, bmax, for a 
region 0.4° < |6| < bmax, 0° < |Z| < l^ax — ^max for the continuum spectrum of a 100 
GeV neutralino annihilating into W^W~ (dot-dashed line, NFW profile) and for the 
gamma-line emission for the benchmark "model I" in the inert Higgs doublet scenario 
(solid lines). From left to right, the solid lines refer to the Moore, NFW and Kravtsov 
profiles, see Table I. From Ref. [59]. 

feeling of the quantitative difference with respect to the GC case treated in the previous 
section, let us note that for our usual benchmark model {rax = 100 GeV, (av) = 1 pb 
with bb final state), a "naive" count statistics above 1 GeV (including all backgrounds) 
would lead to a S/N ~ 6 for the 2° x 2° window considered in Sec. 4, while for the 
region considered above and femax = /max = 25° one would find almost 40000 events from 

§ Even for heavy dark matter particles of TeV mass, the continuum photon spectrum peaks at or below 
30 GeV, explaining in which sense the line emission considered here is "extreme" . 
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DM vs. about seven millions background events, resulting in a S/N > 14. Of course, 
this analysis docs not take into account angular and energy cuts that (as shown for the 
GC case above) do improve the diagnostic power, as well as other possible systcmatics. 
Yet, the importance of this diffuse signal cannot be underestimated (and it is actually 
confirmed by other analyses, see Sec. 4.2 in Ref. [6]), especially in the case no signal is 
revealed from the GC, which may be due to a inner halo profile more cored than NFW. 
The extended signal has indeed a milder dependence on the profile. Also, from Fig. 6 
it is clear that — at least for the simple S/N estimator — it is only important to adopt 
the right angular cut within a factor of two or so in order not to degrade the sensitivity 
by more than~ 20%. But "blindly" focusing on a too narrow window (degree scale) 
around the GC might be overly penalizing; especially for relatively cored profiles, the 
loss in sensitivity may reach a factor of two or more. 

Besides improving the prospects for detection, other advantages in focusing on an 
extended region around the Galactic Center include: 

• enabling an empirical determination of the dark matter profile slope outside of the 
region dominated by the gravitational potential of the central supermassive black 
hole. 

• obtaining independent evidence that an "excess" signal with respect to backgrounds 
is the product of DM annihilation, by comparing the emission from many angular 
regions. In contrast, if the spectrum is found to vary with location, it is most likely 
the product of astrophysical backgrounds. 

• constraining the quantity of dark matter substructure in the halo, by observing the 
angular distribution of the emission. The shape of the smooth, unresolved inner 
profile [22] could be studied, in addition to an anisotropy/multipole analysis [60, 
61, 62]. This, in turn, may have important implications for DM cosmology; see the 
review by T. Bringmann in this issue. 

6. Conclusions 

The most challenging task for indirect dark matter searches is not to detect a few events, 
but to confidently identify those events as the products of dark matter annihilations. In 
particular, any signal must be separated from astrophysical backgrounds if it is to be 
reliably claimed to be a detection of dark matter. This is certainly true in the case of 
gamma-ray telescopes hoping to observe dark matter annihilations in the region of the 
Galactic Center (GC). The discovery of a bright TeV source of astrophysical origin at the 
GC by lACTs has changed the prospects for such searches considerably. One strategy 
in light of this is to focus on different targets with lesser astrophysical background (such 
as dwarf galaxies). Another is to search for a signal in the "noisy" GC region, by taking 
advantage of the pecuhar spectral and angular properties expected from dark matter 
annihilation products. 

In this paper, we have reviewed the latter strategy, reporting on recent studies of 
the GC and inner Galaxy. Given the characteristics of this search, the LAT instrument 
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on board of the Fermi satellite is better suited than existing lACTs. Fermi-LAT 
instrument will detect a number of astrophysical sources in the region of the sky around 
the GC, including the point sources previously identified by HESS and EGRET, and 
perhaps others. A diffuse gamma-ray background will also likely be present. Although 
predictions of Fermi's sensitivity are unavoidably limited by our incomplete knowledge of 
these backgrounds, we have shown that the spectral and angular differences between the 
signal and backgrounds should be distinctive enough to allow one to separate signal from 
background over a significant region of the parameter space, at least for a sufficiently 
cusped dark matter profile (NFW-hke or steeper). 

In the optimistic case where dark matter annihilation products are identified by 
Fermi, then it may also be possible to measure or constrain the properties of the dark 
matter, including its mass, annihilation cross section, and spatial distribution. It is 
unlikely that Fermi will determine the WIMP's mass with high precision, however. For 
example, for the case of a 100 GeV WIMP with an annihilation cross section of 3 x 10^^^ 
cm^/s, distributed with an NFW halo profile, the mass could be determined to he within 
approximately 50-300 GeV. In the same benchmark model, the inner slope of the dark 
matter halo profile could be determined to ~ 10% precision. The combination of several 
indirect detection channels will be crucial to both confirm such a detection, and to best 
constrain the WIMP's properties. On the other hand, it is not excluded that Fermi will 
lead to a radical revision of the present gamma-ray picture of the GC, revealing a more 
complicated zoo of astrophysical accelerators than envisaged in the present estimates. 
In the case where cither the DM signal from the GC is too low or the the background 
is too large/ complex, a DM discovery in gamma rays is still possible by looking at the 
emission from an extended region in the inner halo with Fermi, or from other dark 
matter substructures with both Fermi-LAT and lACTs. In particular, the morphology 
and the spectral properties of the unresolved Galactic background at £■ ^ GeV will be 
useful to optimize the angular and energy-cut templates for searches of the DM emission 
from an annulus of several tens of degrees around the GG. 
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Appendix A. The Case of Decaying Dcirk Matter 

In the case of decaying dark matter, Eq. (1) is modified to 




(A.l) 
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where F is the decay width (inverse hfetime) and the spectrum now refers to the 
photons generated in the decay process. Unhke with the cross section in the case 
of annihilating dark matter, one does not have any strong theoretical motivation for 
considering any particular lifetime for an unstable DM particle. In any case, arguments 
have been put forward justifying the typical range of the lifetimes needed for significant 
signatures in astrophysics with ~TeV mass particles and GUT-scale physics mediating 
the process (in analogy with the expected proton decay in GUTs), see e.g. [63]. From the 
phenomenological point of view, there are a couple of points worth mentioning regarding 
decaying DM candidates: 

lo ^ i I 




Figure Al. The angular profile of the gamma-ray signal as function of the angle, 0, 
to the center of the galaxy for a NFW halo distribution for decaying DM (solid red 
line), compared to the case of self-annihilating DM (dashed blue line). Both signals 
have been normalized to their values at the galactic poles, 9 — ±90°. The central cusp 
is regularized by assuming in both cases an angular resolution of 0.1°. 



I. The DM distribution and the role of substructures in particular is of little 
importance in determining the level of the signal. 

II. The angular distribution of the gamma-ray signal is very distinctive, and much 
flatter than the corresponding annihilation signal, as illustrated for a NFW profile 
in Fig. (Al). 

Should gamma rays be detected from DM, a comparison between the emission in 
the inner Galaxy and the emission at high latitude would immediately reveal the nature 
of the particle physics process (annihilation or decay) responsible for the emission [64]. 
Notice that this information is very difficult to extract with other cosmic ray probes. 
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Appendix B. A Comment on the Factorization Assumption in Eq. (1) 

Although the factorization between the particle physics term and astrophysical term in 
Eq. (1) is a useful approximation and valid in most practical cases, there are exceptions. 
More correctly, one should write 

AAf 1 r [ r 1 

^^{E^,fl) - -T^iE,)- 2- / / a{Vr)VrU{Vr)d^Vr p\i,n)di, (B.l) 

where Vj. is the relative velocity between the two particles (with Vr = |vr.|) and u{vr) its 
distribution function (not necessarily isotropic), whose integral over d^v^ is normalized 
to unity. The factorization assumed in Eq. (1) holds only if the integral in square 
brackets — which is nothing but (av) — is independent of position. A sufficient condition 
for this is that o'(fr) Vr is velocity-independent. In general, the integral depends on 
the kinematical structure of the halo via the position-dependent velocity dispersion, 
anisotropy, etc. Then, both the astrophysical distribution of the DM and the particle 
physics contribute in determining the angular shape of the signal. 

One case in which the factorization is not valid can be found when the WIMP 
annihilations mainly through a P-wave process, such that a{vr)vr oc [65]. This is 
a largely academic case, however, since whenever P-wave annihilation is dominant the 
gamma-ray signal is expected to be suppressed. More interesting is the case in which 
non-perturbative processes lead to large "Sommerfeld enhancements" to the annihilation 
cross section at low velocities [66]. This effect can be thought of as the distortion of 
the wave-function due to a relatively long-range attraction between the WIMPs. In this 
case, a further steepening of the signal towards the Galactic Center is expected, which 
in turn should ease the detection of gamma rays from the inner Galaxy [67]. 
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